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Abstract 

We study the neutrino sector in a minimal SU (3)lxU (1).y model, in which its mass is generated 
at one-loop level with the charged lepton mass, and hence there exists a strong correlation between 
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I. INTRODUCTION 


The flavor problem can be one of the biggest challenging issues to be resolved, since the 
standard model (SM) does not have any theoretical sources (especially) why the number 
of family for each fermion sector (quarks and leptons) is three. One of the reasonable 
interpretations is to extend the gauge sector SU(2)l x U(1)y to SU(3)l x U(1)x, so called 
3-3-1 model, in which the origin of three family is coming from the number of SU(3) color 
of quarks (that has three) due to the gauge anomaly cancellation [l|, 2], Because of larger 
gauge group comparing to the SM one, there are several variations of models extending the 
Higgs sector 3-7} and revisited models to reanalyze with current experimental data |8H29]. 

On the other hand, explaining the current neutrino oscillation data and dark matter 
(DM) candidate might be done by physics beyond SM. It is known that in recent radiative 
seesaw models they can be only explained simultaneously but also correlated each other. 
It means that neutrinos do not directly interact with the SM Higgs field but with a DM 
candidate. As a result, minuscule neutrino masses can be naturally realized. This is because 


a vast literature has recently arisen along thought of this subject 


30 


|4l|. 


In this paper, we combine the 3-3-1 model and radiative seesaw model based on a minimal 
model in Ref. [5|, in which we do not impose any additional discrete symmetry. Then we 
can generate the neutrino mass at one-loop level. Moreover since our model has a strong 
correlation between the charged lepton sector and the neutrino sector due to the same origin 
of these masses, it may be worth analyzing the neutrino oscillation data as well as lepton 
flavor violating processes and so on. 

This paper is organized as follows. In Sec. II, we show our model including Higgs potential. 
In Sec. Ill, we analyze lepton sector and show how to correlate the charged-lepton masses 
and neutrino masses. Then we also show to compute lepton flavor violating processes and 
muon anomalous magnetic moment. In Sec. IV, we perform parameter scan to identify 
allowed parameter regions. We conclude in Sec. V. 
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II. MODEL SETUP 


We discuss a possibility of a one-loop induced radiative seesaw model in the context of 
3-3-1 model in ( 5 ]. The particle contents are shown in Tab. |U We introduce a gauge triplet 
fermion L l = ( u L ,eL,e c R ) with t/(l)x=0. For new bosons, we introduce SU(3)l sextet 
scalars S with U(l)x=0, SU(3)l triplet scalars (r),p,x) with U(l)x = (0, 1 , — 1), respec¬ 
tively. The renormalizable Lagrangian for Lepton Yukawa sector, and the scalar potential 
under these assignments are given by 


C Y = y tl Y L L i(L L ) C jril€ ijk + ye 2 TT[L L S(L L ) c ] + h.c. (II.l) 

i,j,k= 1—3 

V = m 2 |p| 2 + m 2 p \p\ 2 + m\ |x| 2 + m|Tr[|F| 2 ] 

+ ^iM 4 + ^2|p | 4 + A 3 |x | 4 + A 4 [Tr[|5'|]] 4 + A 5 Tr[ A| 4 ] + A 6 |p| 2 |p|' + A 7 |p| 2 |x| 2 + A s |x| \P 
+ Ag|p| 2 Tr[ S\ 2 ] + A 10 |p T| 2 ] + An|x| 2 Tr[|S' 2 ] + A 12 | p^pY + A i 3 |x t? ?| 2 + ^h\p*x\ 2 

+ (mi Yj ^ViPjXk + h.c.) + ( m 2 p T S ] \ + h.c.) + ( m^rf Sp + h.c.) 

ij,k= 1—3 

1-3 1-3 1-3 

+ (777.4 Y Y elJkelmn SiiSjmSkn + h.c.) + (/ 5 ( ? ? t p)(h t x) + h.c.) + (/ 6 Y ^ k SiiPjXkVt + h.c.) 

i,j,k l,m,n hjj* 

1-3 1-3 

+ f 7 Tr[S*Sri*ri T ] + f 8 Tr[S'Sp*p T ] + f 9 Tr[S'S X *X T ] + (fio Y Y ^ lmnS u S jmVkVn + h.c.), 

l,j,k / , 772,71 

( 11 . 2 ) 
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Scalar Fields 
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TABLE I: Contents of lepton and scalar fields and their charge assignment under SU(3)l x U(1)x, 
where the index of the generation are abbreviated. 


3 















where is an anti-symmetric 3 by 3 matrix and yt 2 is a symmetric one. Here the scalar 
helds can be parameterized as 
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(IL4) 


Here v = ^jv £ + + v 2 p = 246 GeV is the vacuum expectation value (vev), and we have 


fl. 


The 


assumed that the vev of cr° is zero, which is unlike the original model discussed in 
scale of v x , which breaks SU(3)l symmetry, is assumed to be of 0(1) TeV, while the other 
three vevs (v a ,v p ,v n ) are assumed to be 0(1-100) GeV. 


III. LEPTON SECTOR 


A. Charged lepton sector 


The charged lepton masses can be generated by the terms ye 1 JV ■ fc=1 _ 3 
and ye 2 Tr[L,LS(L l) c \ after the electroweak symmetry breaking, as can be seen in Eq. (III.II) . 
The resulting mass matrix can be written and diagonalized by 

■*' dia 9 )i(Ve R )ib(e R )b , 


(e L ) a (me)ab(e R ) h = (e L ) a 


(Vh)v v (ye 2 )v a 


L V2 


V2 


(e R ) b = (e L )a(V^ L )ai(m 


ab 


E 


(yh) v v + (y^) v o 


(l/<iK + (Vh) v o 


= (.V^Um^(Vuh, 


L V2 V2 l V2 V2 ibk 
where is the ant-synmretric matrix and p/ 2 is the symmetric matrix. 


(mi) 

(III.2) 


\m 


diag i 


(\m e \, |m M |, \m T \) = (0.5 MeV, 105.6 MeV, 1777 MeV). Each of V e L and V eR is the uni¬ 
tary matrix to diagonalize the charged leptons. From Eq. (IIII.1D and the properties of 
anti-symmetric and of symmetric p/ 2 , the charged lepton Yukawa couplings p/ 1(2) can 
explicitly be rewritten in terms of the charged lepton mass matrix as 


(iJil )ab 


(mi) 


ab 


(mi) T 


ab 


Vzi 


(Vli) ab 


(mj) ab - (m e )j b 

V^Vo- 


(III-3) 
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FIG. 1: An example of one-loop diagrams contributing to neutrino mass. 


Therefore, these Yukawa couplings can be taken as the output parameters, once the charged 
lepton masses and their mixings are fixed as inputs in the following numerical analysis. 


B. Neutrino sector 


Through the following Yukawa interactions, 

C D y h (V L e c L -e L v c L )ih + y h {-u L eR + e c R u c L )rjf 

+ yi 2 {v L e c L + e L v c L )hi + y f „ 2 {p h e R + e R u c L )ht, (III.4) 

the active neutrino mass matrix m v is induced at one-loop level (see Fig. |T] for an example 
diagram), which is given by 
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(III.5) 


Note here that the mass insertion approximation has been used, that is, (M + ) 2 = (M + ) diag + 
(M+) 2 ff _ diag with (M+) 2 iag > (M+) 2 ff _ diag = Sm+ 2 f ,, (/*, f' = rhj), where M + is the singly 
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charged boson mass matrix, and 





m 2 v x v p 

V a 


m 2 + 


rn i v x v p 

v v 


m 2 + 


Sm mm = v^ m 3 v v + (f5V P v x )/2, <5m£ 2 fta = 3\[2m^v a , 


Al3^, 


6m mhi = (4/10 + h) v a v r\l (2v / 2), Sm^ 2 hi = m 3 v v + {f s v p v x )/{2y/2), 
5m S 2 = m 3 ^ + (f 6 v P v x )/(2V2), Sm+ 2 h2 = (4/ 10 + h)v a v v /(2y/2). 


(III.6) 


It suggests mi, m 2 » m 3 , m 4 and / 9 , A 13 » / 5 , / 6 , A, /io- 

Neutrino mass matrix can be diagonalized by a unitary matrix t/j, as m u = U u m^ mg U^, 
and hence the MNS matrix Umns is defined as 


Uum = V',U„. 


(HI-7) 


C. Lepton Flavor Violations 


Here we discuss the lepton flavor violating processes in our model, which can be induced 
through the mixings among singly charged bosons in Eq. (1111.41) . The branching ratio of the 
lepton flavor violating decays are given by 

127T 2 

BR(«, -> tjl) = 2r<~2 [M + l^w.1 2 ] . (III.8) 

rn l ijp 

where 


A-Lji r/T/j 

-4 Rji ITlj 


Gi‘\r]2, h 2 ) + G^ l+V \r]2, h 2 ) 
Gr (^7i, hi) + G ( x + \rji,hi) 


+ m 


+ rrii 


G^irjijhi) + G^ v+m \r)i, hi) 
Gr\t] 2 : h 2 ) + Gft l+V \r]2, h 2 ) 


(III.9) 
(III.10) 


GlIr^Vu hi) ~ 

GfwirkM) ~ 


Sm mhi\.(ylyhY + (yjy^y} 

12(4tt ) 2 m 2 l m 2 i 

5m nh*\(yiiy\y + (yhylj} 

12(47r) 2 m 2 2 m 2 2 
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(iv-\-vi) 

L{R) 


(V2,h 2 
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{yiivly (ye 2 yl 
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V2 
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(111.11) 

(111. 12 ) 

(111.13) 


with (y\„yej = VeRy\ a yi b V} Rl and (yt a y\y = V eL y la y\v} L) with (a, b) = 1,2. The upper 
indices of G, respectively, represent a difference among the singly charged bosons exchanging 
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processes. Therefore, (i) is the r]f and hf mixing process, (ii) is the rj^ and h\ mixing one, 
(Hi) is the 77 ^ no mixing one, (■ iv ) is the 77 .^ no mixing one, (u) is the h\ no mixing one, and 
(vi) is the h\ no mixing one. 

There are box types of LFV precesses i~ —» through y il2 in general, however we 

expect such processes are negligibly small when our mass insertion approximation is reliable. 


The penguin types of these L 


one-loop induced processes 


42 ]. 


'V processes are also negligible, comparing to the G —> ij'y 


D. Muon Anomalous Magnetic Moment 


The muon anomalous magnetic moment, so-called the muon g — 2, has been measured at 
Brookhaven National Laboratory. The current average of the experimental results is given 
by Q 


a® xp = 11659208.0(6.3) x 10 


-10 


(III.14) 


It has been known that there is a discrepancy from the SM prediction by 3 . 2 cr 44] to 


4.1 <7 


45j: 


A^V = C p - a T = (29.0 ± 9.0 to 33.5 ± 8.2) x 10 


-10 


The contribution by the singly charged bosons is evaluated as 


m. 

A an =-- 

m 2 


2 „ 


71 , hi) + Gy ( 77 !, hi) + G [ l ; ( 772 , h 2 ) + G R ; ( 772 , h 2 


kb 


ii), 


■»(**) / 


+2 Gf i+v) 


( 771 , hi) + 2 G ( ' l + \rj 2 ,h 2 ) 


*=3=V 


(III.15) 


(III.16) 


IV. NUMERICAL ANALYSIS 


For simplicity, we assume the hierarchical neutrino mass spectrum and set the lightest 
eigenvalue to be zero. Then we parametrize the neutrino mass matrix as 

M™ p = Umns diag(7n I/1 , m V2 , m U3 ) NS , 

with the (real) standard form of the MNS matrix, 


GmNS = 


1 0 0 


c 13 9 S 13 


C 12 S 12 0 

0 c 23 S 23 


0 1 0 


—S 12 C 12 0 

0 — s 23 C 23 


—S 13 0 c 13 


0 0 1 


(IV. 1) 


(IV.2) 
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with Sij = sin dij and Qj = cos 9 tJ . Then, one finds a relation between M .and A4f, xp given 
by 

MT P = Vj L Mt h V eL , (IV.3) 


where we also define V e u R ) as the standard parametrization as analogy to the MNS matrix, 


14 


eL(R) = 


b C eL(R) 23 S eL(R) 23 
b ~ S eL{R )23 C eL(R )23 


C eL(R) 13 b S eL (R)i 3 


0 1 

:, eL(_R)13 


r eL(R) 

13 


C eL(R) 12 S eL{R) 12 b 
— S eL(i?)12 C eL(i?)12 b 


(IV.4) 


wit h s eL(K) .. = sin 9 eL{R)ij and c eL{R)ij = cos 0 eL{R)ij . 

We perform parameter scan to reproduce the following neutrino oscillation data at 95% 


confidence level 


46] in the next subsection; 


0.2911 < si 2 < 0.3161, 0.5262 < s 2 23 < 0.5485, 0.0223 < s 2 13 < 0.0246, (IV.5) 

| ml 3 -ml 2 1 = (2.44 ± 0.06) x 10" 3 eV 2 , m 2 2 - m 2 Ui = (7.53 ± 0.18) x 10" 5 eV 2 . 

We also impose the current experimental bounds on the lepton flavor violating processes H 
48|: 

BR (fi —> ey) < 5.7 x 10^ 13 , BR(r —> /xy) < 4.4 x 10~ 8 , BR(r —» ey) < 3.3 x 10~ 8 . 

(IV.6) 


A. Normal ordering 

In case of the normal ordering, (m Vl ,m P2 ,m U3 ) = (0, m U2 ,m U3 ), the input parameters vary 
in the following ranges, 

1 GeV < [v v ,v a ) < 100 GeV, 0.1 GeV 2 < (Sm+ 2 hi , Sm+ 2 hi ) < 10 GeV 2 , 

100 GeV < m v + 2 < 1000 GeV, -1 < (s eL{R)ij , s eL(R)lj ) < 1. (IV.7) 

The mass parameters m h + , Sm^ m , and Sm^ 2 h2 can be written in terms of these input pa¬ 
rameters from Eqs. (1111.51) and flIV.ll) . And their mass scales are found to be m h + ~ 0(100) 
GeV, and 5m R 2 h2 ~ <5m+ 2 2 ~ 0(1) GeV 2 . Totally we have twelve input parameters shown 
in Eq. (1IV.7[) to reproduce the three lepton masses (3 charged lepton masses and 2 neutrino 

















mass differences) and three mixings (3 MNS mixings). Once these input parameters are de¬ 
termined, all the physical values are uniquely fixed as discussed in the previous subsections. 
Then we can compare the observable or upper constraints from each of the experimental 
result such as neutrino oscillation data and LFVs. 

Here we show some representative figures of our results which simultaneously satisfy 
the neutrino oscillation data and the constraints from the LFV processes. Here we have 
examined 10 7 sampling points to search for our allowed parameters. The left panel of Fig. [2] 
shows the allowed points in terms of s eL12 and s eL23 to simultaneously satisfy the neutrino 
oscillation data and the LFV constraints, except the red points which predict too large LFV 
rates. We see that that LFV constraints are not so stringent. The left panel of Fig. [3] 
shows the allowed points to satisfy the current LFV bounds in terms of m + , along with the 
future reach in Mu2e experiments 


49] at around BR(/i —> ey) = (2.5 — 6) x 10 17 (horizontal 


lines). For the allowed point, we have also calculated the contribution to the muon g — 2 from 
charged bosons. The left panel of Fig. [4] shows the predictions in terms of m + . We have 

Vl 

found that the contribution is negative and cannot reconcile the discrepancy between the 


experimental result and the SM prediction. See [ 9 ] for the contribution from extra gauged 
bosons, which is found to be positive. 


B. Inverted ordering 


In case of the Inverted ordering, (m ul ,m U2 , m U3 ) = (m ul ,m U2 , 0), we have to fine-tune 
our free parameters to satisfy the neutrino oscillation data, since the neutrino mass matrix 
is nealy diagonal. Hence, the number of solution is very limited compared to the normal 
ordering case. However, the LFV constraints would be automatically satisfied, once we find 
solutions. Such a situation is often seen in supersymmetric models (see, for example) 50]. 
The input parameters vary in the following ranges, 


33 GeV < v v < 37 GeV, 44 GeV < v a < 48 GeV, 0.1 GeV 2 < (Sm^L , ^ 1 GeV ^ 

100 GeV < m + < 1000 GeV, —1 < s e Li 2 < 0, 0 < s e L 23 < 0.5, —1 < s e Li 3 < —0.7, 

' 1 1,2 

— 0.4 < s e Ri2 < —0.1, 0 < s e R 23 < 0.2, 0.4 < s e i?i 3 < 0.7. (IV. 8 ) 

The mass parameters anc ^ ^ m hih 2 can written in terms of those above 

parameters, and their mass scales are, respectively, ~ 0(100) GeV, and 8m^ h2 ~ 
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FIG. 2: Parameter scan in terms of s e _Li 2 and s eL 23 to satisfy the LFV constraints and neutrino 
oscillation data. Here the left panel is for the normal ordering case, where the red points do not 
satisfy the LFV constraints. The right panel is for the inverted ordering case. 


^ m viV 2 ~ (9(1) GeV 2 . Then we show some representative figures of our results which simul¬ 
taneously satisfy the neutrino oscillation data and the constraints from the LFV processes. 
Here we have examined 2 x 10 7 sampling points to search for our allowed parameters. The 
right panel of Fig. [2] shows the allowed points in terms of s eLl2 and s eL23 to simultaneously 
satisfy the neutrino oscillation data and the LFV constraints. We see that that LFV con¬ 
straints do not affect to our solutions as expected. The right panel of Fig. [3] shows the 
allowed points to satisfy the current LFV bounds in terms of m + , along with the future 

Vl 

reach in Mu2e experiments 49] at around BR(/i —> ey) = (2.5 — 6) x 10 -1 ' (horizontal 
lines). Comparing to the normal case, the resulting BR(yU —y ey) is much smaller and more 
difficult to be detected even with the future experiments. For the allowed points, we have 
also calculated the contribution to the muon g — 2 from charged bosons. The right panel of 
Fig. H] shows the predictions in terms of m r] +- One finds that the result is similar to the case 
of normal ordering. 


V. CONCLUSIONS 

We have proposed a radiative seesaw model with a SU{3)c x SU(3)l x U(1)x gauge sym¬ 
metry, in which the neutrino mass is induced through one-loop radiative corrections with the 
charged lepton mass. As a result, there is a strong correlation between the charged lepton 
and neutrino masses, and it is nontrivial if the current neutrino oscillation data are repro- 
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denote the future reach by Mu2e experiments 


49j ] at around BR(^t —»■ ey) = (2.5 — 6) x 10 17 . The 


left panel is for the normal ordering case, while the right panel is for the inverted ordering case. 
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FIG. 4: Muon anomalous magnetic moment in terms of m r) + to satisfy the LFV constraints and 
neutrino oscillation data, where the red points do not satisfy the LFV constraints. One can see 
that all the absolute values of data is less than ICG 10 . Since the observed deviation from the SM 
prediction is positive and of 0(1CF 9 ), the other contributions such as extra gauged bosons could 


be expected 


9(]. Here the left panel is for the normal ordering case, while the right panel is for the 


inverted ordering case. 


duced. In the model, the LFV processes are also induced via one-loop quantum corrections. 
We have performed general parameter scan for the normal and inverted mass ordering cases, 
and found that a large portion of parameter space can simultaneously satisfy the current 
neutrino oscillation data and the constraints on the LFV processes. The parameter region 
we have found can be partly tested by future Mu2e experiments. We have also calculated 
contributions to the muon anomalous magnetic moment from charged scalar particles in our 
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model, and found that the contributions are not significant. 
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